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We consider analytically the evolution of a three-dimensional wave packet generated by 
an impulsive source in a mixing layer. The base flow is assumed to be spreading due 
to viscous diffusion. The analysis is restricted to small disturbances (linearized theory). 
A suitable high-frequency ansatz is used to describe the packet; the key elements of 
this description are a complex phase and a wave action density. It is found that the 
product of this density and an infinitesimal material volume convecting at the local 
group velocity is not conserved: there is a continuous interaction between the base flow 
and the wave action. This interaction is determined by suitable mode-weighted 
averages of the second and fourth derivatives of the base-flow velocity profile. 
Although there is some tendency for the dominant wavenumber in the packet to shift 
from the most unstable value toward the neutral value, this shift is quite moderate. In 
practice, wave packets do not become locally neutral in a diverging base flow (as do 
instability modes), therefore, they are expected to grow more suddenly than pure 
instability modes and do not develop critical layers. The group velocity is complex; the 
full significance of this is realized by analytically continuing the equations for the 
phase and wave action into a complex domain. The implications of this analytic 
continuation are discussed vis-a-vis the secondary instabilities of the packet: very 
small-scale perturbations on the phase can grow very rapidly initially, but saturate later 
because most of the energy in these, perturbations is convected away by the group 
velocity. This remark, as well as the one regarding critical layers, has consequences for 
the nonlinear theories. 
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1 . Introduction 

An unstable shear flow, when excited by an external disturbance that is compact in 
space and impulsive in time, is capable of dispersing this concentrated disturbance into 
a wave train (or wave packet) of finite extent, which is convected by the underlying 
(unstable) base flow. The classical work is due to Gaster (1975; Gaster & Grant 1975). 
The term dispersion is used here in a generalized sense: not only do the various Fourier 
components of the initial disturbance travel at different phase speeds (resulting in a 
broadening of the spatial extent of the disturbance at subsequent instants of time), but 
they also possess different growth rates. Therefore, any concentrated disturbance in the 
flow evolves into a wave packet owing to the 'slippage' of the various Fourier 
components with respect to each other and to the changes in the relative magnitudes of 
these components, because of their different growth rates. 

It is fair to say that our theoretical understanding of how such a wave packet 
evolves in a homogeneous parallel base flow, under the small disturbance approximation 
(i.e. linearized analysis), is quite complete (Criminale & Kovasznay 1962; Tam 1967; 
Gaster 1975; Balsa 1988, 1989a). This is because we may use Fourier-Laplace 

transforms in (x - z) space and time (see figure 1) to solve the initial value problem and 
then evaluate the corresponding inverse transforms, for large values of time, by 
standard asymptotic methods. This yields a concise mathematical description of a wave 
packet, including its receptivity with respect to the disturbance-producing impulse 
(Balsa 1989a). 

On the other hand, when the underlying base flow is inhomogeneous, either in 
space or time, Fourier-Laplace methods are no longer applicable. If the inhomogeneity 
is weak in the sense that the base flow changes slightly over one wavelength of the 
primary oscillation in the packet, a suitable high-frequency analysis may be used (Balsa 
1989b). Before we continue our discussion of this analysis, it is necessary to recall that 
our wave packet has a modal structure (in the ^-direction) and the above-mentioned 
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Fourier components are associated with ( x . z)-space; the convection and dispersion of 
the packet also occur in this space. Therefore, any weak inhomogeneity of the base 
flow is understood to be in (x, z)-space and possibly in time, and the dependence of the 
base flow on the cross-space variable, y, is arbitrary (see figure 1). Finally, we 
emphasize that our attention is focused on wave packets generated by spatially compact 
and temporally impulsive disturbances; it is possible to generate wave packets (or, 
precisely, modulated wave trains) of much longer spatial extent by slowly varying the 
amplitude or frequency of a spatial instability mode that is excited by a wave-maker. 
We plan to examine the properties of these modulated wave trains in a forthcoming 
paper. 

The evolution of spatial instability modes in slightly inhomogeneous base flows has 
been examined by a number of authors; perhaps the most relevant work for our 
purposes is that of Crighton & Gaster (1976). These investigators found that as the 
mode propagates downstream in a slowly diverging jet flow, the instability becomes 
weaker and, as a result, the amplitude of the mode grows at a diminishing rate. In 
fact, there is a streamwise location at which the spatial instability mode is locally 
neutral and, at this location, the modal amplitude levels off. In other words, the 
instability wave saturates because of flow divergence. Although these observations are 
helpful in the interpretation of certain experimental data (e.g. Gaster, Kit & Wygnanski 
1985), Goldstein & Leib (1988) recognized a very different significance, which provided 
the framework for the rational composite description of the nonlinear evolution of an 
instability mode. 

Goldstein & Leib (1988; see also Hultgren 1992) argued that the evolution of an 
instability mode in a diverging base flow merely represented the disturbed flow in a 
suitable outer (i.e. upstream) region. This outer region slowly drives the most unstable 
mode (generated at some location, say, x ■» -oo) toward its neutral state; just prior to the 
attainment of this state, the disturbance develops a distinct critical layer (Maslowe 



- 4 - 


1981). The flow in this layer is dominated by nonlinear effects, hence the outer region 
provides the 'initial' conditions for this nonlinear flow analyzed by Goldstein & Leib 
(1988) using asymptotic methods. 

In a broad sense, our original objectives in embarking on this research were 
twofold: first, the generalization of the ideas of Crighton and Gaster in describing the 
behavior of three-dimensional instability wave packets evolving on top of slightly 
inhomogeneous base flows and, second, the application of critical layer concepts, in the 
spirit of M. E. Goldstein and his colleagues, toward a description of the nonlinear 
phenomena in these packets. The present paper deals with the first of these objectives 
in considerable detail, focusing on physical and mathematical issues relating to the 
possible development of critical layers, the role of the complex group velocity and wave 
action, and the secondary instabilities of the wave packet owing to the presence of an 
inhomogeneous base flow. Several novel effects arise in this study — these have 
important implications for the linear and nonlinear development of wave packets. 

There are many approaches to studying the evolution of rapidly varying 
disturbances (i.e. the wave packet) on a slowly varying system (i.e. the base flow). If 
the disturbance is periodic, the 'method of averaging' yields the equations on the long 
scales; when this procedure is applied to the perturbation Lagrangian, a very powerful 
method results for the description of the disturbance on the slow scales (Whitham 1974; 
Hayes 1970). Although this method has been used for instability waves, it is only 
applicable when the growth rates are asymptotically small (Landahl 1982). 

For instability waves with finite (but numerically small) growth rates, a high- 
frequency ansatz, in terms of complex amplitude and phase, is the proper way to 
proceed (Itoh 1981; Balsa 1989b). One major difference between the analyses in these 
two references is that the latter introduces the concept of a complex modal wave action, 
^t, which is the product of the square of the amplitude and suitable mode-shapes 
integrated over cross-space. The importance of the wave action is clear from the 
classical theories for conservative wave systems. 
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The present paper continues from where Balsa (1989b) left off. The outline is as 
follows: In §2 the problem is formulated, the base-flow and disturbance equations are 

briefly discussed, and the equations for the phase and the wave action are provided for 
reference so that the reader does not have to consult Balsa (1989b). except for details. 
Section 3 deals with the description of a self-similar base flow and the replacement of a 
spatially diverging base flow by a temporally spreading flow. We believe that this 
replacement, which offers a number of technical advantages leading to closed-form 
solutions, is a very good approximation. A more complete justification is provided 
later. In §§4 and 5, results are given for the phase and the wave action. The phase, 
which obeys a canonical dispersion relation because of a Squire-like transformation, 
involves an integral of the local thickness of the mixing layer; this integral accounts for 
history (or age) effects as the wave packet travels downstream, encountering sections of 
a mixing layer whose thickness is slowly increasing. In §5 we provide a simple 
expression for the complex group velocity, show that its divergence is real, and recast 
the equation for the wave action into a form, (266), that is most suitable for the 
discussion that follows in §6. The important physical and mathematical results are 
discussed in several sub-sections of §6; here attention is called to analytic continuation 
(near the end of §6.2) and to the secondary instabilities in §6.3. 

2. Formulation of the problem 

In order to provide some feel for the key equations that will be used in subsequent 
sections, we rely heavily on the general literature on the evolution of instability waves 
in slightly inhomogeneous base flows (see Balsa 1989b and references cited therein) and 
on classical kinematic wave theory (Whitham 1974) and its variants for modal waves 
(Hayes 1970). Of course, with respect to the latter remark, our formalism is necessarily 
very different from what is used in the classical theories because the concept of an 
'average Lagrangian' cannot be used for instability waves. Generally, these waves are 
not periodic in space (or time), except when their growth rates are vanishingly small. 
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2.1. General remarks 

For the unperturbed base flow, consider an incompressible free-shear flow (such as 
a laminar mixing layer) whose velocity and pressure are denoted by U B and p B as 
const., respectively. At moderately large Reynolds numbers, this flow is approximately 
"parallel" in the sense that 

U B - [1 + 0( € 2)]U(j>. I. t) + e[l + <?(e2)]F(y. r)e 2 (la) 

p B - const. + 0(e 2 ) , e « 1 (16) 

where (e] , e 2 , e 3 ) are the unit vectors along the (x, y, z ) directions, respectively; U • e 2 
- 0 (in other words, eV « 1 is the ^-component of U B at lowest order). We now 
proceed to introduce the terminology and notation; the geometry is illustrated in figure I . 

Assume that all variables are nondimensional. Let L ref and t/ ref stand for 
reference quantities that are used to normalize all lengths and velocities, respectively; 
we choose L ni to be the characteristic thickness and £/ ref to be the characteristic speed 
of the mixing layer. Thus, the Cartesian coordinate x - (x, y, z) is nondimensionalized 
by L ref . and time, denoted by t, is normalized by the transit time (Z, ref /£/ ref ). As 
customary in incompressible flows, pressures are measured in terms of p£/Jr ef . where p - 
const, is the density of the fluid. 

The slow spatio-temporal variables used in (la) are defined by 

£ = e(xe, + ze 3 ) , t = ct (lc) 

where e = (/?e)~* is a small parameter. Here. Re - const, denotes the Reynolds number 
of the base flow: Re - £ ref £/ ref /i; » 1 {v - kinematic viscosity of fluid). Note that this 
Reynolds number is based on a length scale associated with the thickness of the mixing 
layer rather than on a typical streamwise length scale. For this reason, the expansion of 
the base flow proceeds in powers of e - (/?e)~* « 1 rather than in those of (/te)~'/ 2 . 
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It is well known that, under the stated assumptions, the horizontal component of the 
fluid velocity at lowest order, U, obeys the boundary layer equations 


DU a 2 u 

Dl dy^ 

where the time differentiation, following the base flow, is defined by 


(2a) 


JD 

Dt 


i* 


u R ■ S- 


d_ 

dx 


while the physical law of mass conservation requires 



(2b) 


(2c) 


The slow variation of the base flow in the streamwise direction is contained in the fact 
that (U, V) depend on (x, z, t) through the (slow) variables £ and t only. Thus, DU/Df 
is actually of <9(t); therefore, a nontrivial balance is attained in (2a) at lowest order. In 
other words, the slow streamwise evolution of the base flow, arising from its small 
acceleration (positive or negative), is caused by viscous diffusion. 

Let us now assume that this base flow is perturbed by an arbitrary disturbance 
whose (perturbation) velocity and pressure fields are denoted by u - u(x, t) and p * 
p(x, t), respectively. The relevant linearized equations for these quantities are 


Momentum: 


Du 

Dr 


+ u 



Continuity: 



0 


^ - + cV2 u 
ox 


(3 a) 
(3b) 


where V 2 denotes the Laplacian operator in the variables (x. y, z). Our interest is in 
the case where the perturbations are the instability waves of the base flow; these waves 
may be considered inviscid when the Reynolds number of the base flow is larger than 
about 500 (Betchov & Szewczyk 1963) and the base-flow profile possesses an inflection 
point. Therefore, we ignore the last term in (3a) and assume that, in principle, U B is a 
known function of its arguments. 
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The following remarks should reassure the reader that the neglect of the viscous 
terms in (3a) does not vitiate any of the major conclusions of this paper. First note that 
viscous effects show up in (2a) and (3a); these, respectively, produce a viscously 
spreading mixing layer and a (locally) small correction to the inviscid instability modes 
of this layer. However, over the long (or global) streamwise scale of interest, namely x 
= O(Re) - <9(r 1 ), both of these viscous effects formally produce an order-one correction 
to the disturbance. 

On the other hand, it is well known from the classical theories of wave motion that 
slight inhomogeneities in a base flow can produce important and subtle effects (e.g. a 
complex group velocity in our case) on the wave-like disturbances sustained by this 
flow. It is precisely some of these new physical effects that we wish to uncover in this 
paper. 

To put it more quantitatively, suppose we replace e in (3a) by (Te), where (for the 
purposes of the asymptotics) T - <9(1). Now we have essentially tagged the viscous 
correction to the instability modes via the parameter T (which, of course, ultimately 
must be set to unity), and separated this correction from the effects of a spreading base 
flow. The latter effects show up as variable coefficients in (3a). When this is done, it 
is found that the coefficient, H, in (10a) contains an extra term arising from the viscous 
term in (3a) (see Balsa 1994; Hultgren 1992). Hence, H contains additively the effects of 
base-flow inhomogeneity and those of viscosity. Therefore, it should be remembered 
that the remarks we later direct at H apply only to the effects of a nonuniform base 
flow. 

To summarize, our objective is to provide a family of inviscid solutions for (3) 
when the base flow varies slowly in the streamwise and spanwise coordinates (x, z), and 
time t , owing to the viscous diffusion of this flow. 


1 II 
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2.2. Fluid displacement 

The so-called primitive variables for the incompressible equations of fluid 
mechanics are the velocity and pressure. These variables are convenient to use in a 
large class of small-disturbance problems, especially when the base flow is a uniform 
stream. On the other hand, when the base flow is dependent on space or time, even if 
this dependence is slow in some sense, the use of these primitive variables is not 
convenient for wave-like problems. This is because, often, there is a subtle energy 
exchange between the base flow and the wave. 

Classically, it is well known that for waves that do not exhibit an exponential 
growth or decay, a most powerful technique for capturing this energy exchange is via a 
variational principle and Whitham's (1974) average Lagrangian, which arises from this 
principle. This suggests that the proper dependent variables to use for instability waves 
traveling on an inhomogeneous base flow are those in which the equations of motions 
may be expressed in terms of a simple variational principle. These dependent variables 
are the fluid particle displacement and pressure; for a disturbance about a known base 
flow, the perturbation displacement is clearly the most useful. In addition, it is simplest 
to consider this displacement as a function of x and / rather than that of a Lagrangian 
(i.e. particle) label and time. This leads to a hybrid Eulerian-Lagrangian description for 
the instability wave, which is especially simple for linearized problems (Andrews & 
McIntyre 1978). 

Therefore, in order to provide the disturbance equations in their purest form, we 
replace the disturbance velocity u by a displacement variable a - a(x, t) such that 


u 




(4) 


A physical interpretation for a « a(x, /) is this: Consider a fluid particle of fixed 
identity that occupies the point x (at time /) in the unperturbed base flow. The position 
of this same particle in the perturbed flow (at time t) is defined to be (x + a). In other 
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words, the particle that would have been at x in the base flow is actually in a slightly 
different position (namely, x + a) in the disturbed flow. 

Before we provide an asymptotic solution of (3) using a high-frequency ansatz, it is 
useful to introduce some terminology. First, we think of the ^-direction as the cross- 
space (in which our disturbance has a mode-like structure) and of the |-space as 
propagation space (in which our disturbance has a wave-like structure). Second, the 
entire effect of base-flow spreading is succinctly contained in F (Balsa 1989b), where 



Finally, note that fF lies in propagation space and represents the (nondimensional) force 
acting on a unit volume of fluid in the boundary-layer approximation. 

2.3. High-frequency ansatz 

The basic variables that we use to characterize a disturbance are the displacement 
and perturbation pressure, (or, p). We expand these (see Balsa 1989b for details) 
according to the usual high-frequency ansatz 

a - a(x, t) - (a(°) + faO) + ...) exp (i^/e) + cc (6) 

where i - and cc stands for the complex conjugate of the term immediately 

preceding it. A similar expansion holds for p\ <t> - 0({, r) is the complex phase of the 
disturbance and ab) - aO) (_y, £, r) for / — 0, 1. We emphasize that <j> is complex, 
although usually we refer to it simply as the phase. 

The physical picture is that the base flow evolves on a 'long or slow' spatio- 
temporal scale that is proportional to the Reynolds number while the local behavior of 
the disturbance is determined by 'short or fast’ spatio-temporal scales characterized by 
the shear layer thickness and associated convection times. On the other hand, the global 
evolution of the disturbance is dictated by the slow scales; hence the relevance of a 
high-frequency ansatz. 
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After substituting (4,6) into the linearized disturbance equations (3), and requiring a 
non-trivial solution for (-)(°), we find the evolution equation for the phase 

- U - £2(V0. & T) (la) 

where £2 ■ £2(k, £. r) is, in principle, a known function of its arguments expressing the 
complex local frequency. a> - £2 - -d<f>/d r, as a function of the complex local wave 
vector, k - V<f>. We call £2 the dispersion relation, and much of stability theory deals 
with the determination of £2 for a variety of flows (Drazin & Reid 1981). More 
generally, a dispersion relation is of the form 0(d<f>/d r, V<p, £, r) « 0; in using (7a), we 
have assumed that 0-0 may be solved explicitly for d<f>/dr. For an incompressible 
mixing layer, there is one such solution that represents the unstable mode of interest. 

Note that V is the gradient operator in propagation space (3/3x was used for this 
operator in physical space) defined by 

v.«, A + «, | (74) 

where ? - (£. ?) - £e ( + fe 3 . Both the frequency and wave vector, (w, k), are local in 
the sense that they depend on £ and r due to the viscous spreading of the base flow. 

Furthermore, having satisfied the dispersion relation via (la), the lowest-order 
solutions (ae< 0 ), p^>) must be proportional to an instability mode that can be supported by 
the base flow. If p m - p m (y, {, r, k) is the pressure mode at wave vector k and co 0 is 
the Doppler-shifted frequency defined by 

u 0 -u 0 (y, l t)--^ + U- V*J (8a) 

then the lowest-order solution is given by 

p(°) - /!({. T)p m (y, i, t . V<f>) m 

where A « -4(£, r) is the slowly varying complex amplitude of the disturbance. 
Similarly. 
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a (0) - /!({. TXQe 2 -i W) . W • e 2 - 0 (9a) 

where 

Q - GO'. t) - -L (96 ) 

W - W(y. fc r) - - ^ /> m (9c) 


and p m on the right-hand sides of (96, c) is to be evaluated at k - V0. 

Two remarks are in order: First p m - p m (y, {. r, k) satisfies a Rayleigh-like 
instability equation at wave vector k; the coefficients of this equation depend 
parametrically on | and r owing to the inhomogeneity of the base flow. Second, the 
separation of the displacement aW (or more precisely the displacement mode) into 
components in cross- and propagation spaces (in terms of Q and W) is done for 
convenience (see (10b)). 

In order to solve for the (•)(*) perturbations, a solvability condition must be 
satisfied. This condition determines the evolution of the amplitude A(i, t) or, more 
properly, the evolution equation of the complex wave action density, jC - Mi, r), on 
the long scales 

+ V • (G M) + HyC = 0 (10a) 

where 

Jt - Mi, r) - A 2 w 0 (C 2 + W • W)d y (106) 

J-OO 

and 


H 



e fr • w dy 


(Q 2 + W • W)d y 


(10c) 


- 13 - 


Recall that &" is defined by (5) and note that H accounts for the viscously 
spreading base flow interacting with the inviscid wave packet. It is precisely this 
factor that vanishes in classical kinematic wave theory— thereby (10a) reduces to a 
conservation equation for the wave action density ji. The complex group velocity, G, 
is defined by 


G 


afflk, s, r) 

dk 


k 


V0 


( 11 ) 


so that ultimately both G and H (both complex) depend on £ and t only. Thus (10a) is 
a partial differential equation for jC in propagation space and time; this equation, 
together with that for the phase, (7a), provides the lowest-order uniformly valid 
solutions for the disturbances on the long spatio-temporal scales. In a linear theory, (7a) 
and (10a) are effectively uncoupled. 

The rest of this paper is devoted to obtaining solutions to these equations under 
reasonable physical and mathematical assumptions, and to examining the secondary 
instabilities of these solutions. We remind the reader that a detailed derivation of the 
previously mentioned evolution equations for 0 and jC may be found in Balsa (1989b). 


3. Temporally spreading base flow 

Although the mathematical problems for the base flow and the disturbance may be 
formulated quite concisely and with total generality (as in §2), the solutions to these 
general problems are completely numerical. What we want to accomplish in this study 
is to gain a thorough qualitative understanding of the physics and to generate semi- 
quantitative results for a three-dimensional wave packet evolving in a spreading base 
flow, such as a mixing layer. Therefore, we shall make a number of additional 
approximations. These will enable us to obtain representative and illuminating closed- 


form results. 
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A wave packet is a wave train of finite extent; because of the instability of the 
base flow and the resultant exponential growth of the disturbance, the dominant part of 
a packet consists of a few oscillations only (Gaster 1975; Balsa 1989a). Thus the overall 
length of a wave packet is very much smaller than the long viscous scales on which the 
base flow is evolving, therefore at each instant of time, the spatial spreading or 
divergence of the mixing layer is unimportant. Of importance is that at two instants of 
time, say, T| and r 2 , the wave packet is at successive streamwise locations, say, and 
£2 » at which the thicknesses of the mixing layer may be very different. It is possible to 
simulate this effect by allowing the thickness to depend on time rather than space. 
This approximation is equivalent to a temporally spreading base flow whose 
characteristic thickness is 5 - 5 (r). 

With this approximation at hand, we write 

U - U(y , Tie, , V - 0 (12a) 

so that the base flow is parallel to the x-axis. U obeys the classical diffusion equation 
(to which (2a) reduces under the stated assumptions), for which we can obtain a class of 
closed-form solutions. For example, if the two external streams of the mixing layer are 
brought into contact at some early time, say, r — -1, the vortex sheet between these 
streams quickly diffuses into a mixing layer of finite thickness and U is expressible in 
terms of the error function whose argument is the similarity variable 



where 

h - S(t) - (I + r)i /2 (12c) 

In view of these preliminary remarks, we write 

U « U(t,) (13) 


where U on the right-hand side of (13) is any reasonable function (e.g. error, tanh, etc.) 
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that approximates the velocity profile of a mixing layer and 5(r) is more or less an 
arbitrary monotonically increasing function of time; 5(0) - 1. We include this 
generalization over (12c) in order to estimate the robustness of our physical model, 
although most of our discussion (unless otherwise stated) will invoke (12c) in order to 
represent the effects of molecular diffusion most realistically. However, for 'turbulent' 
diffusion it is much more reasonable to set 5 (t) = 1 + t (Gaster et al. 1985). 

We are now in a position to provide a more complete picture of the generation of 
the wave packet. Suppose our base flow is perturbed at t « 0 (when 5 - 1) by a 
spatially compact disturbance that is impulsive in time. The concentrated disturbance 
imparted to the flow quickly disperses into a wave train of finite extent, which contains 
the exponentially growing part of the disturbance signal, arising from the instability of 
the mixing layer. The entire wave train is convected downstream (at a convection 
speed, say, U c ) while the disturbance at the center of the packet is growing at the 
maximum growth rate associated with temporal instability. The formation and 
evolution of such a wave packet is well understood in the case of a non-spreading flow 
(i.e. 5 = 1); our objective is to describe this evolution in a self-similarly spreading base 
flow (i.e. 5 - S(t)). 

4. Dispersion relation and the phase 

Under the assumption that the base-flow velocity profile depends only on the 
similarity variable i) - yf 5(r), considerable progress can be made in representing the 
dispersion relation and the modes in terms of certain canonical quantities. Let us 
consider a function P « P(r)) (the canonical pressure) that satisfies a Rayleigh-like 
equation 

<l4 ‘> 

with boundary conditions P -* 0 as if ■+ too. Here, k is a given complex number. In 



order to obtain a nontrivial solution to (14a) we must have 


C - C( k) (146) 

and 

P = P(t), k) (14c) 

On the other hand, the actual pressure mode, p m - p m (y, t, k), satisfies the Rayleigh 
instability equation 


[q - wm 2 


d_ 

dy 


1 

[q - wm 2 



= 0 


(15) 


where the wave vector is k « ( k , i) - Ae ( + le 3 and the dispersion relation is 12 - 
S2(k, r). The boundary condition requires that p m -* 0 as y -* ±oo. 

Upon recalling that t ) = y/b(j) and comparing (14) and (15), we find 


P m O', t, k) - P[r,, 6(t) (k 2 + t 2 y/2] 

(16a) 

- C[6(r) ( k 2 + ! 2 )i/*] 

(166) 


where k • k = k 2 + f 2 . In some sense, (16 a,b) may be thought of as a Squire-like 
transformation; more precisely, these equations express the fact that the pressure mode 
and the complex phase speed, Sl/k, depend on the combined entity ( k 2 + i 2 ) 1 ' 2 rather 
than individually on k and !. 

In view of these remarks, it is easy to see that the canonical quantities P and C 
represent the pressure and complex phase speed of a two-dimensional mode (at complex 
streamwise wavenumber k - k) at the initial instant of time (when 5-1). There 
remains to say a few words about the velocity profile, U = f/(i)), and the canonical 
phase speed, C = C(k). We set 

U(v) = U c + m (17a) 


where 
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U, 


Vy *Ml 


(17 6) 


AU = U { - U 2 > 0 


(17c) 


and /(tj) -» ±1 as i) ■* ±°°. Clearly, the speeds in the two external streams of the mixing 
layer are given by £/, (= 1) and U 2 > 0; see figure 1. 

In order to obtain quantitative results, we let 


/(t?) - -{or 


tanh rj 


erffo) 


-OO < T) < CO 


(18 a) 
(186) 


These velocity profiles are plotted in figure 2; also shown in this figure are the two- 
dimensional growth rates of temporal instability (i.e. /cCj(ic); k real) for each of these 
profiles. Here 

C - C(k) - C r (k) + /C,(i c) (19) 


where C R and denote the real and imaginary parts of C, respectively. It is well 
known that in the temporally unstable range of k, C r = U c . 

Both growth-rate curves in figure 2 are quite close to each other and also to a 
parabolic shape with a neutral wavenumber of 2 kq (s 1) and maximum growth rate of y 
(— 0.1) occurring at the wavenumber Kq (st 0.5). We now replace the actual growth-rate 
curves by a single parabolic approximation (roughly the average of the two curves) and 
analytically continue the latter into the complex plane. Thus, for complex k, 

C(> c) - U c - ^ (K - 2k 0 ) (20) 

K o 

where U c , 7, and k 0 are suitable parameters associated with the instabilities of the base 
flow. The evolution equation for the phase becomes, via (7a) and (166). 

+ §f cwt) (0 { 2 + 0 f V 2 ] = 0 


( 21 ) 



where 0^ = d0/3|, 0j = 30/3f, and C, as a function of its argument, is given by (20). 
We shall solve (21) for the phase, 0 = 0 (|, t). 

After substituting (20) into (21), we obtain a first-order nonlinear partial 
differential equation for 0. The solution of this equation, obtained by standard methods 
(Courant & Hilbert 1966. p. 97) is conveniently written in the form 

4> * <Kb r) = - ^ 0 { (0| + 0f) I/2 /(r) (22a) 

K 0 

where 

J - J{r) - [ ^T'Jdr' (22£) 

Jo 

is the integrated thickness of the mixing layer up to the current time. 

In order to make the phase in (22a) an explicit function of £ and r, we must 
express the local wave vector, V0 - (0|, <fy), as a function of these variables. This is 
done in Appendix A. A discussion of the results for the phase, wave vector, etc., will 
be provided in §6. The role of £/ c , 7, and Xq in determining the dynamic similarity of 
the flow is explained in Appendix B. Note that the phase vanishes at the initial instant 
of time, r * 0, when the base flow is disturbed impulsively. 

Some remarks are in order: First, it is possible to obtain C = C(k) numerically, 

tabulate these results on a grid in the complex x-plane, and then use these tabulated 
values (with some form of interpolation) to solve (21) numerically. This we wanted to 
avoid, at least in the present study, in order to develop some analytical feel for the 
behavior of the phase and its secondary instabilities. Second, the approximation of the 
growth-rate curve by a parabola is quite good (even when 7 is finite) and becomes the 
correct asymptotic limit for many flows in the vicinity of the critical value of an 
externally adjustable parameter that separates regions of stability from those of 
instability. Third, for a mixing layer, the maximum temporal growth rate itself is 
proportional to the velocity difference A U - U { - U 2 \ 7 ~ 0.2k o A£/ is a reasonable 
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approximation for both the tanh and erf profiles (see (18) and (17a)). Finally, J/t is the 
average thickness, say, 3 =* 5(t), and the quantity (r5 /J), which appears throughout the 
analysis, may be interpreted as the ratio of the current (or 'local') and the average 
thickness, 5/3 ^ 1. 


5. The amplitude equation 

With the solution for the complex phase out of the way, we now focus our 
attention on the equation for the wave action, (10a). One objective is to derive an 
interesting expression for H, under the assumption that we have a temporally spreading 
and self-similar base flow, as described in §3. 

We begin with the group velocity, G. After expressing the dispersion relation, Q - 
S2(k, t). in terms of the canonical phase speed C - C(k) via (166), employing (20), (11), 
and our results for V0 (see A 3,4), we obtain 



and, from (23 a), 


V • G - 2 8/J 


(23 a) 


(236) 


It is rather easy to establish (23a) for two-dimensional disturbances (in which case only 
the first term on the right-hand side survives); however, for three-dimensional 
disturbances, considerable algebra is needed. 

Note first that the group velocity is generally complex; the imaginary part of the 
streamwise component is proportional to the growth rate, y, while the imaginary part of 
the spanwise component vanishes. Second, for a non-spreading base flow (i.e. 5/3 « 1), 
the group velocity is actually real, 

G - i/r (24) 


This result (i.e. (24)) is well known from the stationary-phase analysis of classical 
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parallel flows (Gaster 1975). Interestingly enough then, wave packets in a spreading 
base flow differ in an essential way from those in non-spreading flows because of their 
complex group velocities; the physical implications of this will be discussed in §6 in 
terms of the secondary instabilities of the packet. Third, the divergence of the group 
velocity (i.e. V • G - 25// - real) is interpreted as the terra responsible for the 
algebraic time-decay of the complex wave action; specifically, momentarily setting 
H = 0 in (10a) we find 

jCJ 2 - const. - adiabatic invariant (25) 

for an observer (formally) moving with the group velocity. In other words (or still with 
H = 0 ), A ~ /'* , a result generalizing the classical stationary-phase analysis that gives 
A - T~ l for a wave packet in a non-spreading base flow. For example, if 5 ~ t 1 / 2 due 
to the molecular diffusion of the base flow, A ~~ /'* ~ t - 3 / 2 ; therefore, the algebraic 
decay of the amplitude is enhanced owing to the effects of flow spreading. This effect 
is even stronger in a 'turbulent' flow where we might set 5 ~ t so that A ~ r -2 . These 
simple scaling rules, especially when it is realized that £ ~ r for a convected wave 
packet in a spatially spreading base flow, make considerable sense physically. 

Based on these opening remarks, we introduce 

0 - (26 a) 

where, in view of (10a), 5S satisfies 

^ + G • VSS + HG8 - 0 (26 b) 

and H is defined by (10c). We now obtain a more explicit representation for H\ it is 
this coefficient that ultimately contains the most important interaction between the base 
flow and the inviscid wave packet. This interaction is completely absent in all classical 
theories of wave action. 


I 1 1 
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First use (16) in the definitions of Q and W, (9b, c) in order to write these equations 
in terms of the canonical pressure mode P\ replace all ^-derivatives and integrals with 
those involving T), use (7a), (8a), and (166) and define 


and 


m<) 


•r 


Oal 


[cm - um* 


P{q, K)P'(v. K)dr) 


(27a) 


D(Kb 


r 

J-oo 


1 


CC(K) - (/(TOP 


{[P'(v. <c)] 2 + K?PHn, K)}drj 


(276) 


F(k) 


m 

m 


(27 c) 


where the prime denotes differentiation with respect to if. Then from (10c), (5), and 
(12a) we find 


H 


F(m<pj + $*/*] 

6 2 (r) 


(28) 


After some integrations by parts in (27 a, 6) and invoking the Rayleigh equation, 
(14a), we arrive at a representation for F(k), namely. 


F(k) 


<cr ,u > 

m 


where (h) denotes a mode-weighted average of h - A( ij), defined by 


(29a) 


< 6 >- 



*)] 2 

[C(k) - (/(*)]< av 


(296) 


Of course, {•) is a function of a single complex variable, say, k. 

Thus, apart from the scale factor S~ 2 , the coefficient of 55 in (266) is expressible as 
the ratio of the weighted averages of the fourth and second derivatives of the base 
velocity profile, U - £/(tj), evaluated at [5(0^ + 0 2 ) 1 / 2 ]. After substituting (A 3,4) for V0 
into (28), we get 
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H = H(i. t) (30) 

In §6 we discuss the behavior of F , both as F = F(k) and as F = /"(£, r), via the 
substitution k = S(tX 0$ + 0j) 1/2 . 

A feature of (29a) is that F(k) is independent of the velocity ratio {U j - U-^HU j + 
U-£). Therefore, H may be determined, once and for all, for each intrinsic velocity 
profile / = /(»?) (see (17), (18), and Appendix B). It is also important to appreciate the 
somewhat misleading form of (28), which arises because our base flow is self-similar; a 
more penetrating interpretation for 8~ 2 is 25’/S. This is valid whenever the base flow is 
spreading due to molecular diffusion (i.e. 5 - (1 + t) 1 / 2 ); here 8' - dS/dr. Thus, H 
vanishes when 8' is zero (i.e. 8 = const.), and for a self-similar profile this takes place 
as 8 ■* oo. These remarks make a great deal of sense physically. 

Finally, observe that when *c equals the neutral wavenumber (say. tc n ss 2 k 0 ) both 
the numerator and denominator of (29a) possess a first-order pole at ij «* 0, and the 
integral in (298) is to be evaluated by letting the contour of the integration pass below 
this pole in the complex tj-plane. In fact, in this special case, the entire contribution to 
the integrals comes from the pole (i.e. the critical layer) because, for real 17, the 
integrands in (29a) are odd functions of rj. For this reason we can evaluate F(n n ) 
analytically. Thus we get a simple result that provides a check on the numerical results 
presented in §6. 

6 . Discussion of results 

6.1. The phase and wave vector 

From (22a) and (A 1,2) it is seen that both the complex phase, 0, and the local 
wave vector, V0, may be written as products of temporal and spatial factors. Introduce 



X 


(31a) 
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Z - -£ (316) 

T 

r-r (3ic) 

and note that 

0 Z) (32) 

where 5 « 5(71 = J(T)/T is the average thickness and as a function of its arguments, 
can be calculated once and for all. We say that a 'fixed' point in the wave packet 
corresponds to X - (X, Z) = const., for example, the point X * Z - 0 is at the center of 
the packet. 

Observe from (32) that for a non-spreading base flow with 5 - 5 - const. - 1 , the 
phase is proportional to T so that the age of a disturbance is directly measured by its 
phase. On the other hand, for a spreading base flow, the phase varies less strongly 
with time because the average thickness of the layer, 5, is monotonically increasing with 
time. For large values of time, we estimate the phase to be proportional to T and 7"° 
- const., respectively, for laminar' (i.e. 8 ~ T l > 2 ) and 'turbulent' (i.e. 8 ~ 71 diffusions of 
the base flow. 

In accordance with the high-frequency ansatz, (6), the imaginary part of the phase, 
01 - lm( 0 ), is a measure of the magnitude of our exponentially growing disturbance on a 
log-scale. In view of our remarks in the preceding paragraph, the disturbance grows 
linearly initially, followed by a milder growth (such as a square root behavior) for 
larger values of time (all on a log-scale, of course). Clearly, the effect of base-flow 
spreading is to inhibit the growth of the disturbance, as expected intuitively, but in a 
base flow dominated by molecular diffusion, the disturbance continues to grow in 
magnitude for all values of time. This is very different from what happens to a single 
instability mode in a spreading base flow; such a mode always becomes 'neutral' for 
some finite value of time (or streamwise location), as demonstrated by Crighton & 
Gaster (1976). We shall provide a convincing explanation for this difference shortly. 
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In figure 3, we show contour plots of the real and imaginary parts of $ as 
functions of ( X , Z). It is important to recall that these plots are independent of time 
and represent the spatial factor of the phase (see (32)). With reference to figure 3(b), 
we see that Im($) - 4>| is negative in the central portion of this figure; thus the wave 
packet is unstable in this region. Furthermore, we may think of the closed curve 
4 > t » 0 as the boundary of the packet and the exponentially large disturbance near the 
middle of the packet diminishes with distance away from its center (at X = 0). Within 
the footprint of the wave packet, the lines of constant phase, 4> R - const., are practically 
parallel (see figure 3a), thereby indicating a wave structure whose wave fronts are 
perpendicular to the streamwise direction (and, in conjunction with figure 3b, a 
structure whose magnitude diminishes in the X and Z directions). These remarks are 
perfectly consistent with Balsa (1989a), in which a surface plot of a wave packet is 
provided in a constant-thickness mixing layer. 

In figure 4 we show the spatial factors associated with the local wave vector V0; 
these are obtained by pulling out the time factor, Kq/S. from V0. In other words, we 
write 

V0 - (0|, 0 f ) - f K(X, Z) (33 a) 

where K - ( K , L) is the desired spatial factor. 

We note from figure 4 that a large region near the center of the packet is 
dominated by a wave vector K, whose real part is roughly (1, 0). Therefore, in this 
same region, 

Re(V0) s ^ e , (33b) 

so that the dominant wavenumber in the packet is continuously shifting to lower values 
with increasing time. This is clear physically: As the thickness of the mixing layer 
increases, the wave packet tries to adjust its length scale such that it stays in step with 
the local length scale of the base flow. Note that initially (at T - 0 with 5-1) the 


I III 
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dominant wavenumber in the packet is the most unstable wavenumber, k 0 , of temporal 
instability. Of course, as the thickness of the mixing layer increases, the local 
dispersion relation is shifting as well; it is possible to visualize this shift on figure 2(b). 
Roughly, with increasing time, both the maximum growth rate and the range of unstable 
wavenumbers diminish— the growth-rate parabola is effectively ’collapsing' toward the 
origin. A most important point to note is that this phenomenon is happening on the 
local or instantaneous thickness. 5 - 8(T), rather than on the average thickness, 5 - 5(T). 
Thus, at any instant of time, the maximum growth rate and the neutral wavenumber are 
given by 

TW - i • ^ (34a,*) 

We are now in a position to explain why the wave packet does not become neutral 
(or does not saturate) as a result of a spreading base flow caused by molecular 
diffusion. From (33*) and (34*), the ratio of the dominant instantaneous wavenumber at 
the center of the packet to the local neutral wavenumber is 


± 

2kq/8 “ 25 


4“ 


oo 


(35) 


Thus, over the largest possible time interval, the do m inant wavenumber of the packet 
shifts from the most unstable value to (roughly) half-way toward the neutral 
wavenumber. Therefore, a wave packet does not become neutral in a spreading base 
flow; in practical terms, a wave packet is likely to grow more violently than a pure 
instability mode in a given base flow. This is exactly what had been observed by 
Gaster (1981). although his qualitative explanation for this observation is very different 
from ours. 

Since the total wavenumber, (<f>^ + 0j)>/», also plays an important role in our 
analysis, its corresponding spatial factor (i.e. (K • K) 1 / 2 ) is presented in figure 5. Very 
roughly, the real and imaginary parts are proportional to (1 - Z^kq/Ij 2 ) and XK 0 /2y, 
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respectively. It is possible to provide a much more accurate, though still approximate, 
representation of these quantities; this we shall not do for the sake of brevity. 


6.2. Complex wave action 

The equation that governs the evolution of the modal wave action is (26 A); in the 
classical theories, SB (of course, real in this case) would be the adiabatic invariant of the 
system. In the present analysis, SB is no longer a conserved quantity because of the 
presence of a nonzero H . Thus 


d SB 
dr 



V 


j & - - H3 


(36 a) 


where d/dr formally expresses time-differentiation for an observer moving with the 
group velocity. Our objective here is to shed some light on the evolution of SB by 
examining the behavior of //, both in wavenumber space and in propagation space. 
Additionally, in the next sub-section, we discuss the implications of a complex group 
velocity vis-a-vis the secondary instability of the packet. 

Since, from (28), 


5 2 (r) 


(366) 


we start with a discussion of F - F(k ). This complex quantity, as a function of 
(complex) x, has been evaluated numerically from its definition, (29a), for both the tanh 
and the erf profiles. The results are presented in figures 6 and 7. 

Note first from figure 5 that the real part of the spatial factor associated with the 
total wavenumber, ($| + 0^) 1/2 , varies from 0 to 1 while its imaginary part varies from 
-1 to 1. Essentially, the complex argument of F in (36 b) ranges over the Cartesian 
product (0, 3x 0 /2) x (-3x^/2, 3xq/2). In this domain, R e(F) - 5 2 R e(H) is negative so that 
the effect of base-flow spreading is to increase the magnitude of SB; this increase is 
algebraic on the slow viscous time scale, r, and its order of magnitude is very roughly 
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t 4 near the center of the packet. Thus, owing to flow divergence, the disturbance is 
capable of extracting additional energy from the base flow via Reynolds stresses 
evolving on the long spatial and temporal scales. Because of the strong similarity of 
figures 6(a) and 6 (b) for the tanh and erf profiles, we believe that our previous 
conclusion is quite general. As mentioned already, at the neutral point, k - k n , F may 
be evaluated analytically; these values are -8 and -6 for the tanh and erf profiles, 
respectively. These theoretical results are entirely consistent with those given in figure 6. 

In figure 7. we show the behavior of Re{F[5(0j + 0f ) ,/2 ]} in propagation space, X, 
for three values of time (r - 0, 2, and 10). These results are obtained by combining 
those in figures 5 and 6. Because of the temporal factor associated with the argument 
of F in the expression above, Re(7 r ) is time dependent; however, this time dependence 
is relatively slow in comparison with the variation of 6 -2 in (36 b). Thus, near the 
center of the packet, S3' 1 dS3/dr a const. /(I + t), where the real value of this constant 
of proportionality is about 4 at early times and is around 5 at large times. This results 
in an algebraic growth of the wave action S3, as discussed previously. 

In order to understand more fully the evolution of the wave action and the role of 
a complex group velocity, we perform a number of manipulations on the governing 
equation, (26 b). First invoke coordinate transformation (31) and introduce a new time, 

FT - log l(T) when dS/dr > 0 (37a) 


so that (26 b) may be thrown into the form 


3S3 

d&- 


+ 



1 ii 
«o , 


9S3 7 9S3 

dX dZ 


— ^ — H S3 - 0 

5(T)-S(D 


(37 b) 


This transformation fails when the base flow is not spreading since (37a) is 
meaningless. Second, separate (37 b) into its real and imaginary parts and then extend 
this 2x2 system of real equations into the complex plane by writing 
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X - X + iY 


(38a) 


3 r (X, Z, 91 + S5 r {X + iY . Z, 5^ « 5, (T, T, Z, 50 


(386) 


and 


0j(*, Z, 50 - &i(* + iT. Z, 50 - B 2 (X. Y, Z, 5") (38c) 


with similar expressions for H - H R + i7/j - (H R , H x ) + (H l , /7 2 ). Our basic 
unknown now is the complex two-vector B - (B t , B 2 )\ it satisfies the Cauchy-Riemann 
equations. Because of this, we can selectively eliminate d/dX in favor of d/idY and 
rearrange (376) as 


dB 

39’ 


+ X 




.2cl v 
« 0 ’ 



z 



m 

5(71 - 5(71 . 


ff i 


n 2 


,-h 2 

.77, 


B 


0 


(39) 


Finally, it is important to realize that Y has nothing to do with the cross-space 
coordinate y, and any physically meaningful solution of (39) must be evaluated at Y - 0. 

The most important difference between (376) and (39) is that the former equation is 
elliptic (with complex characteristics) whereas the latter is a symmetric hyperbolic 
system (actually hermitian) with real characteristics. In fact, if the characteristics of 
(39) are the level surfaces MX. T, Z, 9~) - const., then ¥ satisfies 


[a!" + x dx + Y 3Y + Z az]* * ^ ar “ 0 


(40) 


and thus the parametric representation of the (±) bicharacteristics of (39) are 


9~ « a (41 a) 

X - X 0 e* (416) 

Y - Y 0 e° ± 2l (i _ 

*0 


(41c) 
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Z-Z 0 e* (41 cO 

where <j (» &~) is the parameter and (Xq. Y 0 , Z 0 ) represents a point on the 
bicharacteristic curve at the initial instant of time, r - T « ST ■ 0. In essence, these 
bicharacteristics may be interpreted as the real rays arising from the complex group 
velocity associated with (26b) (or, actually, more precisely with (37b)). 

In order to obtain the solution at an arbitrary point. P - (A", Y - 0, Z), for a certain 
value of time (say, - o), we construct the two 'backward-running' bicharacteristics in 
order to pick up the initial point, ( X 0 , Y 0% Z 0 ), via (41) (see figure 8 for the geometry). 
At this point the initial value of B is available from the analytic continuation of the 
specified initial data in propagation space. Note that even though Y - 0, Yq ? 0. so that 
in order to obtain the physical solution ai P we must have the analytic continuation of 
the initial data at hand. We shall implicitly assume that this continuation is always 
possible. 

Along each of the (±) bicharacteristics, (39) reduces to a first-order, linear, ordinary 
differential equation for the corresponding characteristic variable, (2?, ± iB 2 ). Each of 
these equations may be solved by elementary methods involving the quadratures of H { 
and H 2 . The general qualitative nature of this solution has been anticipated earlier in 
this sub-section. 

It is important to realize that, because of the complex group velocity, there are 
actually two rays (or bicharacteristics) passing through each point in propagation space. 
The projection of these rays onto propagation space is a radial line; the two rays are 
indistinguishable in this space. 


6.3. Secondary instabilities 

Because the group velocity is complex, (37 b) is a real 2x2 elliptic system for 65 ■ 
(£5 r , 65,). For such a system, the initial value problem atT-T-^-Ois 'ill posed' 
in the sense of Hadamard, although for analytic initial data, a properly posed solution 



- 30 - 


exists, at least locally in the immediate neighborhood of the initial hypersurface. We 
interpret the global ill-posedness as the secondary instability of the wave packet. 

In the preceding sub-section we have transformed this ill-posed elliptic system into 
a well-posed hyperbolic system. This transformation is somewhat misleading because 
the analytic continuation of the initial data is, in itself, an unstable process (more about 
this later); specifically a typical rapid, but perhaps spurious, oscillation of wavenumber 
riK 0 (n » 1) along the T-axis is converted into an exponentially large value when 
analytically continued into the complex plane, beyond the strip | Y\ - 0(1 /wcg). The 
two crucial heuristic questions are these: (a) How far do we need to continue the 
initial data into the complex plane in order to describe the evolution of the complex 
wave action for all values of time? (b) What is the role of the group velocity in 
advecting spurious disturbances away from the footprint of the packet, which, after all, 
is confined to a fairly small region in ( X , Z) space (see figure 3)? We will be able to 
answer both of these questions quite satisfactorily by considering the secondary 
instabilities of the packet. 

In order to accomplish this, we perturb the phase by a small quantity. 0, so that 

Total Phase -*•0 + 0 (42a) 

where <f> is given by (22a). Since the total phase satisfies dispersion relation (7a), it is a 
simple matter to show by linearization that the perturbation phase, 0, satisfies 

+ G • V0 - 0 (426) 

This equation is exactly the same as (36a) with H ~ 0, so that the previous 
discussion on analytic continuation is applicable to (426). For example, if \(if\X. Y, Z) 
U - 1. 2) represent, respectively, the analytic continuation of the initial data, \fif\X. Z) 
(if - R - real; if « I * imaginary parts), then by the method of characteristics outlined in 
§6.2, the solution for the analytic continuation of 0 R -*• 0, is 
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1 (X, Y, z, <D = 




A similar equation holds for \ji 2 (i.e. ■* \p 2 )\ these equations are valid for all values of 

time. 

The formal similarity between (43) and the classical D'Alembert solution for the 
simple wave equation is absolutely striking. Apparently the initial data split into two 
halves; each of these propagates along one of the bicharacteristics, and these two halves 
add (or subtract) at subsequent instants of time to form the required solution. 

We are now ready to answer the two questions raised previously. For large values 
of time, ■* oo, the physically valid solution at a fixed (X, Z) is 



(44a) 


(446) 


Since is a logarithmic time scale, (see (37a)), the asymptotic result represented in (44) 
is actually attained only algebraically fast on the physical time scales r or T . 
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The physical interpretation of (44) is this: An arbitrary initial disturbance in the 

phase is advected away from the footprint of the wave packet by the real part of the 
group velocity. This advection is radially outward in propagation space so that any 
disturbance that is not at X - (X, Z) - 0 will be converted out of the packet. 
Therefore, the solution for large values of time will be dominated by the initial 
disturbance near X - 0; this disturbance cannot be adverted away since the real part of 
the group velocity vanishes at the origin. More precisely, it is the analytic continuation 
of this initial data into the interval | Y\ £ 2y/ic 0 that actually matters. Therefore, a 
spurious rapid oscillation of the phase at a characteristic wavenumber, okq (n » 1), will 
attain the value exp(2ny) for very large values of time. 

This result is interesting. It says that perturbations of the phase will saturate; 
therefore, at large values of time, the total phase will be dominated by <fi. This is exact 
in the linear theory. In other words, as ■* oo, ^ -*■ (- complex const.) so that the 

complex phase, <f>, jumps by fixed amounts as a result of the asymptotic response to the 
perturbations in the phase. 

On the other hand, it is easy to see from (43) that for small (or moderate) values of 
time, the above-mentioned rapid oscillations in the phase will grow as exp(2 nyr). This 
observation, combined with the fact that the actual problem is nonlinear, strongly 
suggests that once these small-scale, rapidly growing (though only initially) disturbances 
are triggered in the flow, they will lead to the secondary instabilities of the wave 
packet. 

It is also possible to study the evolution of small disturbances riding on top of the 
wave action, 55, by writing 


Total Wave Action ■* 55 + /5 (45) 

where 5 is the perturbation. After substituting (45) and (42a) into (26b) and linearizing, 
we arrive at an equation for /S. This equation is very similar to (26 b); essentially H is 
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replaced by derivatives of F and O - 3J2/3k with respect to the wave vector, because 
of the perturbations in the phase. It is beyond the scope of this paper to study the 
solutions of this equation; we plan to examine this equation in a forthcoming paper in 
which we will study the secondary instabilities of an instability mode (rather than a 
wave packet). 

7. Conclusions 

It is well known that an instability mode evolving in a spreading base flow attains 
a neutral state at some streamwise location. This is not the case for a wave packet: 
The dominant wavenumber in a packet is such that the local growth rate associated 
with it is in the vicinity of the maximum growth rate of the local dispersion relation. 
Although there is some shift toward the neutral wavenumber (and away from the most 
unstable one), this shift is far from complete even as r (or Q •* oo. For this reason, a 
wave packet is expected to grow more violently than an instability mode in a given 
spreading base flow. 

The wave action, £9, is no longer a conserved quantity; because of the interaction 
between the base flow and the disturbance, S3 grows algebraically on the slow time 
scale, t. The group velocity is complex, and the equation for the wave action may be 
converted into a symmetric hyperbolic system via analytic continuation. We use this 
hyperbolic system to study the behavior of the perturbed phase. It is found that small- 
scale disturbances to the phase grow very rapidly initially and then saturate at later 
values of time. We interpret this initial growth as the secondary instability of the 
packet. 
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stability. 


Appendix A. Dependence of wave vector on space and time 

We obtain a formal solution to (21) using the characteristic differential equations 
belonging to (21). Since the wave vector. V0, is constant along the complex rays, the 
trajectory of these rays may be solved for quite simply, with the initial condition that 
all rays emanate from the spatially compact source of the disturbance at { - 0. We 
express the equations of the rays in (£, r) space as 


I - ~ ^/Kn 

o 


- + tf ) 1 ' 2 - 






02 ) 1/2 


L_ _ hh. 

nJ/x-l (#f + tfyn 


(A 1) 
(A 2) 


where the wave vector V0 - (0 { . 0 f ) may be thought of as the label that identifies a 

specific ray in propagation space (or, more precisely, in the analytic continuation of this 
space). 

In order to obtain V0 as a function of £ and r. we must invert (A 1,2). Let a and 
b temporarily denote the left-hand sides of (A 1,2). respectively. We first form the 
ratio {a/b\ rearrange this to obtain a quadratic equation for (0 f /0 { ) whose solution, in 
conjunction with (A 1,2), yields 


0. _ [q 2 + 462 + a5] i/i 
* 2v/2 
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S - (a 2 - SPY'* (A 6) 

The outside square root on the right-hand side of (A 5) is chosen such that lm(0) - 
0j < 0 (see (6) and (22a); thus we have an exponential instability), while the other 
fractional powers in (A 3,4) and (A 6) are chosen so that, in the first place, 0$ and 0 j. 
are analytic functions of | - (£, J) and that, second, the spatial derivatives of 0, as 
obtained from (22a), are equal to (A 3,4). Although it is cumbersome to verbally 
express these constraints on the fractional powers, it is simple to program them on a 
computer by examining the alternatives in 'if statements.' 

In essence, (A 3,4) provide V0 as functions of ({, r) and. when these functions are 
substituted into (22a), we obtain 0 - 0(£, t). This is the desired result for the phase. 

Appendix B. Dynamic similarity and scale transformation 

The theory outlined in the main body of this paper contains three free parameters 
(see 20): the convection speed. U c — (t/ j + 4/ 2 )/2; the maximum growth rate associated 
with the canonical dispersion relation, y a 0.2x0(471 - i/ 2 ); and the corresponding 
temporal wavenumber, Kq se 0.5. The numbers just provided for y and kq are 
reasonably accurate for both the tanh and erf profiles (see (17) and (18)). 

What we wish to emphasize here is that the dependence on these three parameters 
can be eliminated completely by the introduction of new space variables. Very briefly, 
from (A 1,2), (31), and (A 3,4), we see that 

V0 - ^ (spatial factor) (B 1) 

where the (spatial factor) in (B 1). essentially K of (33a). depends only on the variables 
Xk 0 h and Ztc^/y. Note that X and Z have physical units of speed, and a most 
important speed that totally governs dynamic similarity is the 'instability speed,' y/K$ ~ 
A U = (4/, - U 2 ). 
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Similar remarks hold for the spatial factor of the phase, (32), and the group 
velocity, (23a), relative to the convection speed 



Once again we see from (B 2) the importance of the instability speed, 7 /kq, in the 
normalization of the group velocity. Clearly, any physical mechanism of 'dispersion' for 
an instability wave packet is necessarily very different from that for a classical packet 
in a conservative system. Because of this scale transformation of X » ( X , Z), all the 
results presented in figures 3-5 are essentially valid for arbitrary £/ c , y, and ic 0 . 
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Figure I. Geometry of problem (not to scale), y - cross space; (x, z) - propagation space. 



Figure 2. 



Figure 3. Contours of constant phase, 4>, in the wave packet (U i - 1, t/ 2 - 0, 7 - 0.1, k 0 - 0.5): (a) <t> R - const. 
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Figure 4. Contours of constant slreamwise, K, and spanwise, L, wavenumbers in the wave packet (£/| 
7 - 0.1, Kq - 0.5): (q) )f R - const.; (A) Ky - const.; (c) Z, R - const.; (d) Ly - const. 
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Continued. 
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Propagation space 
(X. Z) 



Figure 8. Bicharacteristics (rays) in the analytic continuation of propagation space. 
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